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I. INTRODUCTION 



Understanding theoretically the patterns of 
masses and mixings of ordinary fermions is one of 
the long aimed goals in particle physics. Of the 56 
parameters in the Standard Model (SM) Yukawa sec- 
tor (including Majorana neutrinos) only 22 can be 
measured at low energy and just 17 have been deter- 
mined from the experiment. Grand Unified Theories 
(GUTs), by enforcing stringent relations among the 
different particle sectors and by reducing the degen- 
eracy in the parameter space, do provide a powerful 
tool for addressing the multiplicity of matter states 
and the detailed structure of the Yukawa sector. 

Appealing candidates for realistic GUTs are mod- 
els based on the SO(10) gauge group [l|. All the 
known SM fermions plus three right-handed neu- 
trinos fit into three copies of the 16-dimcnsional 
spinorial representation of 5*0(10), thus providing 
a rationale for the SM hypercharge structure. The 
model also provides a natural explanation for the 
sub-eV light neutrino masses via the seesaw mecha- 
nism 

The purpose of this paper is to review the con- 
straints enforced by gauge unification on the in- 
termediate mass scales in the non-supersymmetric 
5O(10) GUTs, a needed preliminary step for as- 
sessing the structure of the multitude of the dif- 
ferent breaking patterns before entering the de- 
tails of a specific model. Eventually, our goal is 
to envisage and examine scenarios potentially rel- 
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evant for the understanding of the low energy mat- 
ter spectrum. In particular those setups that, al- 
beit non-supersymmetric, may exhibit a predictiv- 
ity comparable to that of the minimal supersym- 
metric 50(10), scrutinized at length in the last few 
years [4J. 

The most recent discussion of fermion masses and 
mixings in non-supersymmetric 50(10) GUTs was 
given in Ref . Q . The authors focussed only on renor- 
malizable models (i.e. without the spinorial 16# 
in the Higgs sector) with combinations of 10h and 
126ff or 120ir driving the Yukawa interactions. Par- 
ticular attention is paid to the leptonic sector and 
the mechanism of generation of neutrino masses via 
see-saw. 

The constraints imposed by the absolute neutrino 
mass scale on the position of the B — L threshold, 
together with the proton decay bound on the unifica- 
tion scale Mjj, provide a discriminating tool among 
the many 50(10) scenarios and the corresponding 
breaking patterns. These were studied at length 
in the eighties and early nineties, and detailed sur- 
veys of two- and three-step 50(10) breaking chains 
(one and two intermediate thresholds respectively) 
are found in Refs. @, 0, 0, Q . 

We perform a systematic survey of 50(10) unifi- 
cation with two intermediate stages. In addition to 
updating the analysis to present day data, this reap- 
praisal is motivated by (a) the absence of U{1) mix- 
ing in previous studies, both at one- and two-loops 
in the gauge coupling renormalization, (b) the need 
for additional Higgs multiplets at some intermediate 
stages, and (c) a reassessment of the two-loop beta 
coefficients reported in the literature. 

The outcome of our study is the emergence of 
sizeably different features in some of the breaking 
patterns as compared to the existing results. This 
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allows us to rescue previously excluded scenarios. 
All that before considering the effects of threshold 
corrections [13, EH, EH > that are unambiguously as- 
sessed only when the details of a specific model are 
worked out. 

It is remarkable that the chains corresponding to 
the minimal 50(10) setup with the smallest Higgs 
representations (10#, 45# and 16#, or 126# in the 
renormalizable case) and the smallest number of pa- 
rameters in the Higgs potential, are still viable. The 
complexity of this non-supersymmctric scenario is 
comparable to that of the minimal supersymmetric 
50(10) model, what makes it worth of detailed con- 
sideration. 

In Sect. [II] we set the framework of the analysis. 
Sect. Mil provides a collection of the tools needed for 
a two-loop study of grand unification. The results 
of the numerical study are reported and scrutinized 
in Sect. [iVl Perspectives for further progress are 
discussed in Sect. [Vj Finally, the relevant one- and 
two- loop /3-coefficients are detailed in Appendix [A] 



II. THREE-STEP 50(10) BREAKING 
CHAINS 

The relevant 50(10) -> G2 -> Gl -> SM sym- 
metry breaking chains with two intermediate gauge 
groups G2 and Gl are listed in Table HI Effective 
two-step chains are obtained by identifying two of 
the high-energy scales, paying attention to the possi- 
ble deviations from minimality of the scalar content 
in the remaining intermediate stage (this we shall 
discuss in Sect. HVB]) . 

For the purpose of comparison we follow closely 
the notation of ref. [9|, where P denotes the unbro- 
ken D-parity [l3[ . For each step the Higgs represen- 
tation responsible for the breaking is given. 

The breakdown of the lower intermediate symme- 
try Gl to the SM gauge group is driven either by 
the 16- or 126-dimcnsional Higgs multiplcts 16h or 
126#. An important feature of the scenarios with 
126# is the fact that in such a case a potentially 
realistic 50(10) Yukawa sector can be constructed 
already at the renormalizable level. Together with 
10 h all the effective Dirac Yukawa couplings as well 
as the Majorana mass matrices at the SM level 
emerge from the contractions of the matter bilin- 
ears 16^16^ with 12Qh or with 16# 16#/A, where A 
denotes the scale (above Mjj) at which the effective 
dimension five Yukawa couplings arise. 

The Higgs transforming as 10 under 50(10) may 
carry in general extra quantum numbers of a com- 
plex representation of some additional symmetry 
(a discussion on the implementation of a Peccei- 
Quinn U(1)pq symmetry in this scenario is given 



(^Vi p\ in 




G2 




Gl 


I: 


y 

210 


{2 L 2 R A C } 


> 

45 


{2i2 fl l x 3 c } 


II: 


54 


{2 L 2 fl 4 c P} 


210 


{2 l 2r1 x 3 c P} 


III: 


54 


{2 L 2 fl 4 c P} 


45 


{2 L 2 fl lx3 c } 


IV: 


210 


{2 L 2 R 1 X 3 C P} 


45 


{2 L 2 R l x 3c} 


V: 


210 


|2l2r4 c > 


45 


\2 L \ R A C } 

[_ *■* L-j XT A U J 


VI: 


54 


{2 L 2 fl 4 c P} 


45 


{2 L lfl4c} 


VII: 


54 


{2 L 2 fl 4 c P} 


210 


{2 L 2 fl 4c} 


VIII: 


45 


{2 L 2 fl lx3 c } 


45 


{2 L l fl lx3 c } 


IX: 


210 


{2 L 2 fl lx3 c P} 


45 


{2 L l fl lx3 c } 


X: 


210 


{2 L 2 fl 4c} 


210 


{2 L l fl lx3 c } 


XI: 


54 ' 


{2 L 2 fl 4 c P} 


210 


{2 L l fl lx3 c } 


XII: 


45 


{2 L lfl4c} 


45 1 


{2 L l fl lx3 c } 



TABLE I: Relevant 5O(10) symmetry breaking chains via 
two intermediate gauge groups Gl and G2. For each step the 
representation of the Higgs multiplet (in 50(10) notation) re- 
sponsible for the breaking is given. The breaking to the SM 
group ly2^3 c is obtained via a 16 or 126 Higgs representa- 
tion. The naming and ordering of the gauge groups follows 
the notation of ref. Q. 

in Ref. ||). In this case it is sufficient to consider 
only two complex symmetric matrices Yio and Y126 
at the renormalizable 5O(10) level, namely 

16 F (y 10 10 H + Yi 26 1 26 ff )16 F , (1) 

that govern all the effective Yukawa couplings at 
lower energies. Such scenarios are rather constrained 
and hence their detailed numerical studies are well 
motivated . 

D-parity is a discrete symmetry acting as charge 
conjugation in a left-right symmetric context [13j |. 
and as that it plays the role of a left-right symmetry 
(it enforces for instance equal left and right gauge 
couplings). 50(10) invariance then implies exact 
D-parity (because D belongs to the 50(10) Lie al- 
gebra). D-parity may be spontaneously broken by 
D-odd Pati-Salam (PS) singlets contained in 210 or 
45 Higgs representations. Its breaking can therefore 
be decoupled from the SU(2)r breaking, allowing 
for different left and right gauge couplings. 

The possibility of decoupling the D-parity break- 
ing from the scale of right-handed interactions is 
a cosmologically relevant issue. On the one hand 
baryon asymmetry cannot arise in a left-right sym- 
metric (gL = 9r) universe On the other hand, 
the spontaneous breaking of a discrete symmetry, 
such as D-parity, creates domain walls that, if mas- 
sive enough (i.e. for intermediate mass scales) do 
not disappear, overdosing the universe [lj|. These 
potential problems may be overcome either by con- 
fining D-parity at the GUT scale or by invoking infla- 
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Surviving Higgs multiplets in 5*0(10) subgroups 
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TABLE II: Scalar multiplets contributing to the running of the gauge couplings for a given S*O(10) subgroup according to 
minimal fine tuning. The survival of </> 126 (not required by minimality) is needed by a realistic leptonic mass spectrum, as 
discussed in the text (in the 2^2ijlx3 c and 2l1 r 1x3 c stages only one linear combination of <f> 10 and <f> 12e remains). The 
U(l)x charge is given, up to a factor \/3/2, by (B — L)/2 (the latter is reported in the table). For the naming of the Higgs 
multiplets we follow the notation of Ref. [ij with the addition of <^> 126 . When the D-parity (P) is unbroken the particle content 
must be left-right symmetric. D-parity may be broken via P-odd Pati-Salam singlets in 45h or 210h- 



tion. The latter solution implies that domain walls 
are formed above the reheating temperature, enforc- 
ing a lower bound on the D-parity breaking scale of 
10 12 GeV. Realistic 5O(10) breaking patterns must 
therefore include this constraint. 



A. The extended survival hypothesis 

Throughout all three stages of running we assume 
that the scalar spectrum obeys the so called ex- 
tended survival hypothesis (ESH) [l(| which requires 
that at every stage of the symmetry breaking chain 
only those scalars are present that develop a vacuum 
expectation value (VEV) at the current or the subse- 
quent levels of the spontaneous symmetry breaking. 
ESH is equivalent to the requirement of the mini- 
mal number of fine-tunings to be imposed onto the 
scalar potential [TtJ so that all the symmetry break- 
ing steps are performed at the desired scales. 

On the technical side one should identify all the 
Higgs multiplets needed by the breaking pattern un- 
der consideration and keep them according to the 
gauge symmetry down to the scale of their VEVs. 
This typically pulls down a large number of scalars 
in scenarios where 126jy provides the B — L break- 
down. 

On the other hand, one must take into account 
that the role of 126# is twofold: in addition to trig- 
gering the Gl breaking it plays a relevant role in the 
Yukawa sector (Eq. ijl])) where it provides the nec- 
essary breaking of the down quark - charged lepton 
mass degeneracy. For this to work one needs a rea- 
sonably large admixture of the 126# component in 



the effective electroweak doublets. Since (2,2, l)io 
can mix with (2, 2, 15)j26 onr y below the Pati-Salam 
breaking scale, both fields must be present at the 
Pati-Salam level (otherwise the scalar doublet mass 
matrix does not provide large enough components of 
both these multiplets in the light Higgs fields). 

Note that the same argument applies also to the 
2l1.r4c intermediate stage when one must retain the 
doublet component of 126/f, namely (2, +i, 15)j26", 
in order for it to eventually admix with (2, +1, l)io 
in the light Higgs sector. On the other hand, at 
the 2i2flljsf3 c and 2il/jljjf3 c stages, the (minimal) 
survival of only one combination of the 4> 10 and 126 
scalar doublets (see Table [TTJ) is compatible with the 
Yukawa sector constraints because the degeneracy 
between the quark and lepton spectra has already 
been smeared-out by the Pati-Salam breakdown. 

In summary, potentially realistic renormalizable 
Yukawa textures in settings with well-separated 
5*0(10) and Pati-Salam breaking scales call for an 
additional fine tuning in the Higgs sector. In the 
scenarios with 126r/, the 10w bidoublet (2, 2, l)io, 
included in Refs [a, 0, H, must be paired at 
the 2i2fl4c scale with an extra (2, 2, 15)y26 scalar 
bidoublet (or (2, +1, l)i with (2, +1, 15)^ a t the 
2l1/j4c stage). This can affect the running of the 
gauge couplings in chains I, II, III, V, VI, VII, X, XI 
and XII. 

For the sake of comparison with previous stud- 
ies [1, 0, 1, [1 we shall not include the 126 multi- 
plets in the first part of the analysis. Rather, we 
shall comment on their relevance for gauge unifica- 
tion in Sect. lTVCl 
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III. TWO-LOOP GAUGE 
RENORMALIZATION GROUP EQUATIONS 



In this section we report, in order to fix a consis- 
tent notation, the two-loop renormalization group 
equations (RGEs) for the gauge couplings. We con- 
sider a gauge group of the form f/(l)i® f/(l)jv® 
Gi ® ... ® Gjv'j where Gi are simple groups. 



while the other terms stand for the two-loop cor- 
rections, including that induced by Yukawa interac- 
tions. The latter is accounted for in terms of a factor 

Y*{F P ) = ^-yTr [C 2 (F p )YY^ , (6) 

where the "general" Yukawa coupling 

Y abc i>J> h h c + h.c. (7) 



The non-abelian sector 



Let us focus first on the non-abelian sector corre- 
sponding to Gi <£>...<£> Gjv' and defer the full treat- 
ment of the effects due to the extra U(l) factors to 
section UlI B[ Defining t = log(/i//io) we write 



includes family as well as group indices. The cou- 
pling in Eq. |(7J| is written in terms of four-component 
Weyl spinors ip a ,b and a scalar field h c (be complex 
or real) . The trace includes the sum over all relevant 
fermion and scalar fields. 
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B. The abelian couplings and U(l) mixing 



where p = 1,...,N' is the gauge group label. Ne- 
glecting for the time being the abelian components, 



the /^-functions for the Gi 
read at two-loop level [HI, [1 



P 



-y^ 2 (G p ) 



(4tt) 2 

2k y(f) + g2p 

— i n{Fp) + Wf 

20 



G y g auge couplings 
20, 2J: 



4 1 

-kS 2 (F p ) + -r}S 2 (S p ) 



(4tt) 2 
4C 2 (F P ) 



34 9 
-y(G 2 (G p )) 2 



3 
2 



G 2 (G P ) K S 2 (F P ) 



(3) 



4G 2 (S P ) + -G 2 (G P ) ) rjSziSp) 



KC 2 {F q )S 2 {F p ) + T)C 2 (S q )S 2 (S p ) 



where k = 1, \ for Dirac and Weyl fermions respec- 
tively. Correspondingly, r\ = 1, | for complex and 
real scalar fields. The sum over q ^ p corresponding 
to contributions to f3 p from the other gauge sectors 
labelled by q is understood. Given a fermion Fora 
scalar S field that transforms according to the rep- 
resentation R = Ri (g) ... (g) Rn>, where R p is an irre- 
ducible representation of the group G p of dimension 
d(R p ), the factor S 2 (R P ) is defined by 



S 2 (R P )=T{R P ) 



d(R) 
d{R p ) 



(4) 



where T(R p ) is the Dynkin index of the represen- 
tation R p . The corresponding Casimir eigenvalue is 
then given by 



In order to include the abelian contributions to 
Eq. ((3|) at two loops and the one- and two-loop effects 
of U(l) mixing [22| . let us write the most general 
interaction of N abelian gauge bosons A^ and a set 
of Weyl fermions ipf as 



(8) 



The gauge coupling constants g r b, r,b = 1,...,N, 
couple to the fermionic current J p = ip^^Q^ipf. 
The N x N gauge coupling matrix g r b can be diago- 
nalized by two independent rotations: one acting on 
the U(l) charges Q r j and the other on the gauge bo- 
son fields A^. For a given choice of the charges, g r b 
can be set in a triangular form (g r b — for r > b) by 
the gauge boson rotation. The resulting N(N+ 1) /2 
entries are observable couplings. 

Since F* in the abelian case is itself gauge invari- 
ant, the most general kinetic part of the lagrangian 
reads at the renormalizable level 
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P uF a F bfll/ 



(9) 



where a ^ b and \P, a b\ < 1- A non-orthogonal rota- 
tion of the fields A£ may be performed to set the 
gauge kinetic term in a canonical diagonal form. 
Any further orthogonal rotation of the gauge fields 
will preserve this form. Then, the renormalization 
prescription may be conveniently chosen to main- 
tain at each scale the kinetic terms canonical and 
diagonal on-shell while renormalizing accordingly 
the gauge coupling matrix g r b . Thus, even if at 



C 2 (Rp)d(R p ) = T{R p )d(Gp 



(5) 



where d(G) is the dimension of the group. In Eq. 
the first row represents the one-loop contribution 



1 Alternatively one may work with off-diagonal kinetic terms 
while keeping the gauge interactions diagonal [23j . 



5 



one scale g r b is diagonal, in general non-zero off- 
diagonal entries are generated by renormalization 
effects. One shows [2Jj that in the case the abelian 
gauge couplings are at a given scale diagonal and 
equal (i.e. there is a £7(1) unification), there may ex- 
ist a (scale independent) gauge field basis such that 
the abelian interactions remain to all orders diagonal 
along the RGE trajectory 2 . 

In general, the renormalization of the abelian part 
of the gauge interactions is determined by 



dg r b 
~~dT 



— 9raftab • 



(10) 



where, as a consequence of gauge invariance, 



with Z3 denoting the gauge-boson wave-function 
renormalization matrix. In order to further simplify 
the notation it is convenient to introduce the "re- 
duced" couplings [H| 



9kb = Qk9rb , 

that evolve according to 
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The index k labels the fields (fermions and scalars) 
that carry £7(1) charges. 

In terms of the reduced couplings the /3-function 
that governs the £7(1) running up to two loops is 
given by [Sill H| 
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where repeated indices are summed over, labelling 
fermions (/), scalars (s) and £7(1) gauge groups (c). 
The terms proportional to the quadratic Casimir 
C 2 (R P ) represent the two-loop contributions of the 
non abelian components G q of the gauge group to 
the £7(1) gauge coupling renormalization. 

Correspondingly, using the notation of Eq. (fH?)) . 
an additional two-loop term that represents the 
renormalization of the non abelian gauge couplings 



2 Vanishing of the commutator of the /3-functions and their 
derivatives is needed 12511 . 



induced at two loops by the £7(1) gauge fields is to 
be added to Eq. ((3|), namely 



A(3 P = 
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k gj c S 2 (F p ) + 77 g 2 sc S 2 (S p ) . (15) 



In Eqs. lfT4 |) -([T5 |) . we use the abbreviation / = F p 
and s = S p and, as before, k = 1, i for Dirac and 
Weyl fermions, while 77 — 1, ^ for complex and real 
scalar fields respectively. 



Some notation 



When at most one U(l) factor is present, and 
neglecting the Yukawa contributions, the two- loop 
RGEs can be conveniently written as 



da. 



dt 
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(16) 



where on = gj /Air. The /3-coefficients and bij for 
the relevant 50(10) chains are given in Appendix [XI 

Substituting the one-loop solution for Oj into the 
right-hand side of Eq. (|16p one obtains 



2tt 



t^ 1 
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where luj = ajOtj(0)/(2Tr) and bij — bij/aj . The 
analytic solution in (fTT)) holds at two loops (for 
ojjt < 1) up to higher order effects. A sample of the 
rescaled /3-coefficients b^ is given, for the purpose of 
comparison with previous results, in Appendix |A1 

We shall conveniently write the /3-function in 
Eq. (fT4| . that governs the abelian mixing, as 
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(18) 



where ^ sr include both one- and two-loop contribu- 
tions. Analogously, the non-abelian beta function in 
Eq. ©, including the £7(1) contribution in Eq. (fl"5j) . 
is conveniently written as 



ftp 



(Anf 



7p- 



(19) 



The 7 P functions for the 50(10) breaking chains con- 
sidered in this work are reported in Appendix IA 11 

Finally, the Yukawa term in Eq. ([6]), and corre- 
spondingly in Eq. (|14p . can be written as 



Yi(F p )=yp k Tv Y fe rJ 



(20) 



where are the "standard" 3x3 Yukawa matri- 
ces in the family space labelled by the flavour index 
k. The trace is taken over family indices and k is 
summed over the different Yukawa terms present at 
each stage of 5O(10) breaking. The coefficients y p k 
are given explicitly in Appendix I A 21 
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D. One-loop matching 



The matching conditions between effective the- 
ories in the framework of dimensional regulariza- 
tion have been derived in [U HH . Let us consider 
first a simple gauge group G spontaneously broken 
into subgroups G p . Neglecting terms involving loga- 
rithms of mass ratios which are expected to be sub- 
leading (massive states clustered near the thresh- 
old 3 ), the one- loop matching for the gauge couplings 
can be written as 



12tt 



C 2 {G) 
12tt 



(21) 



Let us turn to the case when several non-abelian 
simple groups G p (and at most one U(l)x) spon- 
taneously break whilst preserving a U(1)y charge. 
The conserved U(l) generator Ty can be written in 
terms of the relevant generators of the various Car- 
tan subalgebras (and of the consistently normalized 
T x ) as 



Ty = PiTi 



(22) 



where Y^Pi = 1> an d i runs over the relevant p (and 
X) indices. The matching condition is then given by 



-i C2\Gi 



12tt 



(23) 



where for i = X, if present, C 2 = 0. 



Consider now the breaking of N copies of U(l) 
gauge factors to a subset of M elements U(l) (with 
M < N). Denoting by T n [n = 1, N) and by T m 
(m = 1, M) their properly normalized generators 
we have 



T — P T 



(24) 



with the orthogonality condition P mn Pm'n = $mm'- 
Let us denote by g na {n,a = 1,...,N) and by g m b 
(to, b — 1, M) the matrices of abelian gauge cou- 
plings above and below the breaking scale respec- 
tively. By writing the abelian gauge boson mass 
matrix in the broken vacuum and by identifying the 
massless states, we find the following matching con- 
dition 



(25) 



Notice that Eq. ([25)1 depends on the chosen basis 
for the 17(1) charges (via P) but it is invariant un- 
der orthogonal rotations of the gauge boson fields 



3 An early discussion of thresholds effects in SO(10) GUT is 
found in pjj|. 



(gO T Og T = gg T ). The massless gauge bosons A^ 
are given in terms of Ai^ by 



4A» — 



(26) 



where m = 1, M and n = 1, N . 

The general case of a gauge group U{\)\ ® ... ® 
U(1)n <8 G\ ® ... ® Gn' spontaneously broken to 
U(l)x ® ... ® U(1) M with M < N + N' is taken 
care of by replacing (g<7 T ) _1 in Eq. (|25]) with the 
block-diagonal (N + N ) x (N + N') matrix 
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thus providing, together with the extended Eq. (|24l 
and Eq. (|25[) . a generalization of Eq. 



IV. NUMERICAL RESULTS 

At one-loop, and in absence of the U(l) mixing, 
the gauge RGEs are not coupled and the unification 
constraints can be studied analytically. When two- 
loop effects are included (or at one-loop more than 
one U(l) factor is present) there is no closed solution 
and one must solve the system of coupled equations, 
matching all stages between the weak and unification 
scales, numerically. On the other hand (when no 
U(l) mixing is there) one may take advantage of 
the analytic formula in Eq. (fT7|) . The latter turns 
out to provide, for the cases here studied, a very 
good approximation to the numerical solution. The 
discrepancies with the numerical integration do not 
generally exceed the 10 -3 level. 

We perform a scan over the relevant breaking 
scales Mu, M 2 and Mi and the value of the grand 
unified coupling a\j and impose the matching with 
the SM gauge couplings at the Mz scale requiring 
a precision at the per mil level. This is achieved by 
minimizing the parameter 



S = 



„th 



(28) 



where denote the experimental values at Mz and 
af 1 are the renormalized couplings obtained from 
unification. 

The input values for the (consistently normalized) 
gauge SM couplings at the scale Mz = 91.19 GeV 
are |2g| 



ai = 0.016946 ±0.000006, 
a 2 = 0.033812 ±0.000021, 
a 3 = 0.1176 ±0.0020, 



(29) 



corresponding to the electroweak scale parameters 



a,. 



127.925 ±0.016. 



sm 2 w = 0.23119 ±0.00014. 



(30) 
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All these data refer to the modified minimally sub- 
tracted (MS) quantities at the Mz scale. 

For 0(1,2 we shall consider only the central values 
while we resort to scanning over the whole 3er do- 
main for as when a stable solution is not found. 

The results, i.e. the positions of the intermediate 
scales Mi , M 2 and Mjj shall be reported in terms of 
decadic logarithms of their values in units of GeV, 
i.e. m = log 10 (Mi/GeV), n 2 = log 10 (Af 2 /GeV), 
njj = log 10 (M(//GeV). In particular, nu, ni are 
given as functions of n\ for each breaking pattern 
and for different approximations in the loop expan- 
sion. Each of the breaking patterns is further sup- 
plemented by the relevant range of the values of ajj. 



A. U(l) R x U(l) x mixing 

The chains VIII to XII require consideration of the 
mixing between the two U(l) factors. While U(1)r 
and U(l)x do emerge with canonical diagonal ki- 
netic terms, being the remnants of the breaking of 
non-abelian groups, the corresponding gauge cou- 
plings are at the onset different in size. In general, 
no scale independent orthogonal rotations of charges 
and gauge fields exist that diagonalize the gauge in- 
teractions to all orders along the RGE trajectories. 
According to the discussion in Sect. IIII1 off-diagonal 
gauge couplings arise at the one-loop level that must 
be accounted for in order to perform the matching 
at the Mi scale with the standard hypercharge. The 
preserved direction in the Q R,X charge space is given 

by 

QY = s ll QR+ {i QX > (31) 

where 

Q R = I 3R and Q x = yj\ {^) ■ (32) 

The matching of the gauge couplings is then ob- 
tained from Eq. (12"5)) 

g? = P fa*)- 1 1* , (33) 

with 



and 




SB Hi- For all other cases, in which only one J7(l) 
factor is present, the matching relations can be read 
off directly from Eq. (JUJ and Eq. (f23f . 



B. Two-loop results (purely gauge) 

The results of the numerical analysis are organized 
as follows: Fig. [1] and Fig. [2] show the values of njj 
and n2 as functions of ri\ for the pure gauge run- 
ning (i.e. no Yukawa interactions), in the 126# and 
16# case respectively. The differences between the 
patterns for the 126# and 16^ setups depend on the 
substantially different scalar content. The shape and 
size of the various contributions (one-loop, with and 
without U(l) mixing, and two-loops) are compared 
in each figure. The dissection of the RGE results 
shown in the figures allows us to compare our re- 
sults with those of Refs. @, 0, B @- 

Table UTT1 shows the two- loop values of a^ 1 in the 
allowed region for n%. The contributions of the ad- 
ditional (j> 126 multiplets, and the Yukawa terms are 
discussed separately in Sect. IIVCI and Sect. IIVD1 
respectively. With the exception of a few singular 
cases detailed therein, these effects turn out to be 
generally subdominant. 

As already mentioned in the introduction, two- 
loop precision in a GUT scenario makes sense once 
(one-loop) thresholds effects are coherently taken 
into account, as their effect may become compara- 
ble if not larger than the two loop itself (the argu- 
ment becomes stronger as the number of intermedi- 
ate scales increases). On the other hand, there is 
no control on the spectrum unless a specific model 
is studied in details. The purpose of this work is 
to set the stage for such a study by reassessing and 
updating the general constraints and patterns that 
5*0(10) grand unification enforces on the spread of 
intermediate scales. 

The one and two-loop /3-coefficients used in the 
present study are reported in Appendix Table HXl 
in the appendix shows the reduced 6^ coefficients for 
those cases where we are at variance with Ref . Q ■ 

One of the largest effects in the comparison with 
Refs. d, 0, H, [I] emerges at one-loop and it is due to 
the implementation of the U(l) gauge mixing when 
U(1)r ® U(l)x appears as an intermediate stage of 
the 50(10) breaking 4 . This affects chains VIII to 
XII, and it exhibits itself in the exact (one-loop) flat- 
ness of n% , njj and ajj as functions of ri\ . 

The rationale for such a behaviour is quite simple. 



When neglecting the off-diagonal terms, Eq. ([551) 
reproduces the matching condition used in Refs. [6|, 



4 The lack of abelian gauge mixing in Ref. Q was first ob- 
served in Ref. |26j |. 
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6 8 10 

(k) Chain Xlla 



FIG. 1: The values of njj (red/upper branches) and ri2 (blue/lower branches) are shown as functions of n\ for the pure gauge 
running in the 126h case. The bold black line bounds the region n\ < ni. From chains la to Vila the short-dashed lines 
represent the result of one-loop running while the solid ones correspond to the two-loop solutions. For chains Villa to Xlla 
the short-dashed lines represent the one-loop results without the U(l)x ® U(1)r mixing, the long-dashed lines account for the 
complete one-loop results, while the solid lines represent the two-loop solutions. The scalar content at each stage corresponds 
to that considered in Ref. §, namely to that reported in Table [IT] without the </> 126 multiplets. For chains I to VII the two-step 
50(10) breaking consistent with minimal fine tuning is recovered in the rt2 — » np limit. No solution is found for chain Xa. 
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n 2 n v n 2 nu n 2 n u 

16.0 r [ 17 [ 




(g) Chain Vllb (h) Chain VHIb (i) Chain IXb 



n 2 nu n n n 2 iiu 

16 r 



15.5 r 

15.0; 
14.5 • 
14.0 
13.5 : 



15 ■ 



14 - 



' "1 13.0 r ' «1 I 

6 8 10 12 4 6 8 10 12 14 4 6 8 10 

(j) Chain Xb (k) Chain Xlb (1) Chain Xllb 




FIG. 2: Same as in Fig. [T]fbr the 16h case. 



When considering the gauge coupling renormaliza- 
tion in the 2l1r1x3 c stage, no effect at one-loop 
appears in the non-abelian /3-functions due to the 
abelian gauge fields. On the other hand, the Higgs 
fields surviving at the 2lIrIx3 c stage, responsible 
for the breaking to Iy^lSc, are (by construction) SM 
singlets. Since the SM one-loop /3-functions are not 



affected by their presence, the solution found for 712, 
nu and ot\j in the n\ = rii case holds for ri\ < as 
well. Only by performing correctly the mixed \r\x 
gauge running and the consistent matching with ly 
one recovers the expected ri\ flatness of the GUT 
solution. 

In this respect, it is interesting to notice that the 
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Chain 




Chain 




la 


[45.5,46.4] 


lb 


[45.7, 44.8] 


Ha 


[43.7, 40.8] 


lib 


[45.3,44.5] 


Ilia 


[45.5, 40.8] 


Illb 


[45.7, 44.5] 


IVa 


[45.5, 43.4] 


IVb 


[45.7, 45.1] 


Va 


[45.4, 44.1] 


Vb 


[44.3, 44.8] 


Via 


[44.1,41.0] 


VIb 


[44.3, 44.2] 


Vila 


[45.4, 41.1] 


Vllb 


[44.8, 44.4] 


Villa 


45.4 


VHIb 


45.6 


IXa 


42.8 


IXb 


44.3 


Xa 




Xb 


44.8 


XIa 


38.7 


Xlb 


41.5 


Xlla 


44.1 


Xllb 


44.3 



TABLE III: Two-loop values of a^ 1 in the allowed region 
for n\. From chains I to VII, a^ 1 is ri\ dependent and its 
range is given in square brackets for the minimum (left) and 
the maximum (right) value of ri\ respectively. For chains VIII 
to XII, a^ 1 depends very weekly on n\ (see the discussion on 
f/(l) mixing in the text). No solution is found for chain Xa. 

absence of £7(1) mixing in Refs. [1, 0, H, Q makes 
the argument for the actual possibility of a light 
(observable) U(1)r gauge boson an "approximate" 
statement (based on the approximate flatness of the 
solution) . 

One expects this feature to break at two-loops. 
The SU(2)l and SU(3) C /3-functions are affected at 
two-loops directly by the abclian gauge bosons via 
Eq. (fl5|) (the Higgs multiplets that are responsible 
for the U (\)r®U(1)x breaking do not enter through 
the Yukawa interactions). The net effect on the non- 
abelian gauge running is related to the difference 
between the contribution of the U(\)r and U(l)x 
gauge bosons and that of the standard hypercharge. 
We checked that such a difference is always a small 
fraction (below 10%) of the typical two-loop contri- 
butions to the SU(2)l and SU(3) C /3-functions. As 
a consequence, the n\ flatness of the GUT solution 
is at a very high accuracy (10 -3 ) preserved at two- 
loops as well, as the inspection of the relevant chains 
in Figs. [TH21 shows. 

Still at one-loop we find a sharp disagreement in 
the n\ range of chain Xlla, with respect to the result 
of Ref. [Ij]. The authors find rii < 5.3, while strictly 
following their procedure and assumptions we find 
n\ < 10.2 (the updated one- and two-loop results 
are given in Fig. [T]i). As we shall see, this differ- 
ence brings chain Xlla back among the potentially 
realistic ones. 

As far as two-loop effects are at stakes, their rel- 
evance is generally related to the length of the run- 
ning involving the largest non-abelian groups. On 
the other hand, there are chains where ni and nu 
have a strong dependence on n\ (we will refer to 
them as to "unstable" chains) and where two-loop 



corrections affect substantially the one-loop results. 
Evident examples of such unstable chains are la, IVa, 
Va, IVb, and Vllb. In particular, in chain Va the 
two- loop effects flip the slopes of ni and njj, that im- 
plies a sharp change in the allowed region for ri\ . It 
is clear that when dealing with these breaking chains 
any statement about their viability should account 
for the details of the thresholds in the given model. 

In chains VIII to XII (where the second interme- 
diate stage is IiXrXx'&c-, two-loop effects are mild 
and exhibit the common behaviour of lowering the 
GUT scale (nu) while raising (with the exception of 
Xb and XIa,b) the largest intermediate scale (712). 
The mildness of two-loop corrections (no more that 
one would a-priori expect) is strictly related to the 
(ni) flatness of the GUT solution discussed before. 

Worth mentioning are the limits ~ njj and 
ri\ ~ ni. While the former is equivalent to ne- 
glecting the first stage G2 and to reducing effec- 
tively the three breaking steps to just two (namely 
5O(10) — > Gl -> SM) with a minimal fine tuning 
in the scalar sector, care must be taken of the lat- 
ter. One may naively expect that the chains with 
the same G2 should exhibit for n\ ~ n<i the same 
numerical behavior (5O(10) -> G2 -> SM), thus 
clustering the chains (I,V,X), (II,III,VI,VII,XI) and 
(IV, IX). On the other hand, one must recall that the 
existence of Gl and its breaking remain encoded in 
the G2 stage through the Higgs scalars that are re- 
sponsible for the G2— >G1 breaking. This is why the 
chains with the same G2 are not in general equiv- 
alent in the n% ~ ni limit. The numerical features 
of the degenerate patterns (with ni ~ nu) can be 
crosschecked among the different chains by direct 
inspection of Figs. HH2] and Table EH 

In any discussion of viability of the various sce- 
narios the main attention is paid to the constraints 
emerging from the proton decay. In non supersym- 
metric GUTs, this process is mediated by baryon 
number violating gauge interactions, inducing at low 
energies a set of effective dimension 6 operators that 
conserve B — L. In the 50(10) scenarios we consider 
here, such gauge bosons are integrated out at the 
unification scale, and therefore proton decay con- 
strains nu from below. The present experimental 
limit T p (p — > e+7r°) > 1.6 x 10 33 years [29[ implies 

(?2pj io 2 (™^i5) > 5 .2, (36) 

that, for a^ 1 = 45 yields nu > 15.4. Taking the 
results in Figs. [T]-[5] and Table IIIII at face value the 
chains Vlab, Xlab, Xllab, Vb and Vllb should be 
excluded from realistic considerations. 

On the other hand, one must recall that once 
a specific model is scrutinized in detail there can 
be larg e threshold corrections in the matching [T(| 
[TTI . \12\ , that can easily move the unification scale 
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by a few orders of magnitude (in both directions). 
In particular, as a consequence of the spontaneous 
breaking of accidental would-be global symmetries 
of the scalar potential, pseudo-Goldstone modes 
(with masses further suppressed with respect to 
the expected threshold range) may appear in the 
scalar spectrum, leading to potentially large RGE 
effects [30] . Therefore, we shall follow a conservative 
approach in interpreting the limits on the intermedi- 
ate scales coming from a simple threshold clustering. 
These limits, albeit useful for a preliminary survey, 
may not be sharply used to exclude marginal but 
otherwise well motivated scenarios. 

Below the scale of the B — L breaking, processes 
that violate separately the barion or the lepton num- 
bers emerge. In particular, AB = 2 effective interac- 
tions give rise to the phenomenon of neutron oscil- 
lations (for a recent review see Ref. [3l|). Present 
bounds on nuclear instability give r^ucl > 10 32 
years, which translates into a bound on the neutron 
oscillation time r„_ ft > 10 s sec. Analogous limits 
come from direct reactor oscillations experiments. 
This sets a lower bound on the scale of AB = 2 non 
supersymmetric (dimension 9) operators that varies 
from 10 to 300 TeV depending on model couplings. 
Thus, neutron-antineutron oscillations probe scales 
far below the unification scale. In a supersymmetric 
context the presence of AB = 2 dimension 7 opera- 
tors softens the dependence on the B — L scale and 
for the present bounds the typical limit goes up to 
about 10 7 GeV. 

Far more reaching in scale sensitivity are the 
AL = 2 neutrino masses emerging from the see-saw 
mechanism. At the B — L breaking scale the A R 26 
(S R 6 ) scalars acquire AL = 2 (AL = 1) vacuum ex- 
pectation values (VEVs) that give a Majorana mass 
to the right-handed neutrinos. Once the latter are 
integrated out, dimension 5 operators of the form 
vIvlHH 7 generate light Majorana neutrino states 
in the low energy theory. 

In the type-I seesaw, the neutrino mass matrix m„ 
is proportional to YnM^Y^v 2 where the largest 
entry in the Yukawa couplings is typically of the or- 
der of the top quark one and Mr ~ Mi . Given a 
neutrino mass above the limit obtained from atmo- 
spheric neutrino oscillations and below the eV, one 
infers a (loose) range 10 12 GeV < Mi < 10 14 GeV. 
It is interesting to note that the lower bound pairs 
with the cosmological limit on the D-parity breaking 
scale (see Sect. HH) . 

In the scalar-triplet induced (type-II) seesaw the 
evidence of the neutrino mass entails a lower bound 
onjhe VEV of the heavy SU(2) L triplet in 126 H (or 
in 16#16#). This translates into an upper bound on 
the mass of the triplet that depends on the structure 
of the relevant Yukawa coupling. If both type-I as 
well as type-II contribute to the light neutrino mass, 
the lower bound on the Mi scale may then be weak- 



ened by the interplay between the two contributions. 
Once again this can be quantitatively assessed only 
when the vacuum of the model is fully investigated. 

Finally, it is worth noting that if the B — L break- 
down is driven by 126# , the elementary triplets cou- 
ple to the Majorana currents at the renormalizable 
level and m„ is directly sensitive to the position of 
the Gl -> SM threshold M\. On the other hand, 
the ni-dependence of m„ is loosened in the b-type of 
chains due to the non-renormalizable nature of the 
relevant effective operator 16f16_f16h16h/A, where 
the effective scale A > Mjj accounts for an extra 
suppression. 

With these considerations at hand, the constraints 
from proton decay and the see-saw neutrino scale 
favor the chains II, III and VII, which all share 
2 L 2 R ^ C P in the first 5O(10) breaking stage @]. 
On the other hand, our results rescue from obliv- 
ion other potentially interesting scenarios that, as 
we shall expand upon shortly, are worth of in depth 
consideration. In all cases, the bounds on the B — L 
scale enforced by the see-saw neutrino mass excludes 
the possibility of observable U(1)r gauge bosons. 



C. The <j) 12b Higgs multiplets 

As mentioned in Sect. Ill Al in order to ensure a 
rich enough Yukawa sector in realistic models there 
may be the need to keep more than one SU(2)l 
Higgs doublet at intermediate scales, albeit at the 
price of an extra fine-tuning. A typical example is 
the case of a relatively low Pati-Salam breaking scale 
where one needs at least a pair of SU (2) l <B> SU [2)r 
bidoublets with different SU(4)c quantum numbers 
to transfer the information about the PS breakdown 
into the matter sector. Such additional Higgs mul- 
tiplets are those labelled by </> 126 in Table ITT1 

Table IIVI shows the effects of including <fi 126 at 
the SU(A)c stages of the relevant breaking chains. 
The two-loop results at the extreme values of the 
intermediate scales, with and without the 4> 12e mul- 
tiplet, are compared. In the latter case the complete 
functional dependence among the scales is given in 
Fig. [1] Degenerate patterns with only one effective 
intermediate stage are easily crosschecked among the 
different chains in Table [TVl 

In most of the cases, the numerical results do 
not exhibit a sizeable dependence on the additional 
(2, 2, 15)^ ( or +i; 15)j26") sca l ar multiplets. The 
reason can be read off Table [X] in Appendix [5] and 
it rests on an accidental approximate coincidence of 
the cj) 126 contributions to the SU(4) C and SU(2) LtR 
one-loop beta coefficients (the same argument ap- 
plies to the 2^1 r4c case). 

Considering for instance the 2l2rAq stage, one 
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Chain 
la 

Ha 

Ilia 

Va 

Via 

Vila 

XIa 

Xlla 



n u 



9.50 
8.00 
10.5 
10.5 
9.50 
9.50 
11.0 
10.1 
11.4 
11.2 
11.3 
10.5 
3.00 
3.00 
3.00 
3.00 



10.0] 
9.50] 
13.7] 
13.7] 
13.7] 
13.7] 
11.4] 
11.2] 
13.7] 
13.7] 
13.7] 
13.7] 
13.7] 
13.7] 
10.8] 
10.5] 



[16.2 
[10.4 
[15.4 
[15.4 
[16.2 
[16.2 
[11.0 
[10.1 
[14.4 
[14.5 
[15.9 
[16.5 
[13.7 
[13.7 
[10.8 
[10.5 



10.0] 
16.2] 
13.7] 
13.7] 
13.7] 
13.7] 
14.4] 
14.5] 
13.7] 
13.7] 
13.7] 
13.7] 
13.7] 
13.7] 
10.8] 
10.5] 



[16.2 
[18.0 
[15.4 
[15.4 
[16.2 
[16.2 
[15.9 
[16.5 
[14.4 
[14.5 
[15.9 
[16.5 
[14.8 
[14.8 
[14.6 
[14.7 



17.0] 
16.2] 
15.1] 
15.1] 
15.1] 
15.1] 
14.4] 
14.5] 
14.9] 
14.9] 
14.9] 
15.0] 
14.8] 
14.8] 
14.6] 
14.7] 



[45.5 
[30.6 
[43.7 
[43.7 
[45.5 
[45.5 
[45.4 
[32.5 
[44.1 
[40.8 
[45.4 
[33.3 
[38.7 
[36.0 
[44.1 
[39.8 



46.4] 
45.5] 
40.8] 
37.6] 
40.8] 
37.6] 
44.1] 
40.8] 
41.0] 
38.1] 
41.1] 
38.1] 
38.7] 
36.0] 
44.1] 
39.8] 



TABLE IV: Impact of the additional multiplet </> 126 (second 
line of each chain) on those chains that contain the gauge 
groups 2r J 2fl4c or 2^1^4(7 as intermediate stages, and whose 
breaking to the SM is obtained via a 126j^ representation. 
The values of ri2, ny and a^ 1 are showed for the minimum 
and maximum values allowed for rii by the two-loop analysis. 
Generally the effects on the intermediate scales are below the 
percent level, with the exception of chains la and Va that are 
most sensitive to variations of the /3-functions. 



obtains Aa4 = 
l 



x 4 x T 2 (15) = f , and Aa 2 = 
3 x 30 x T 2 (2) = 5, that only slightly affects the 
value of a [/ (when the PS scale is low enough), but 
has generally a negligible effect on the intermediate 
scales. 

An exception to this argument is observed in 
chains la and Va that, due to their 7i 2 ,[/(ni) slopes, 
are most sensitive to variations of the /3-coefficients. 
In particular, the inclusion of (f> 126 in the la chain 
flips at two-loops the slopes of n 2 and njj so that 
the limit n 2 = njj (i.e. no G2 stage) is obtained for 
the maximal value of n\ (while the same happens 
for the minimum n\ if there is no 126 ). 

Fig. [3] shows three template cases where the (j> 126 
effects are visible. The highly unstable Chain la 
shows, as noticed earlier, the largest effects. In chain 
Va the effects of 4> 126 are moderate. Chain VII is the 
only "stable" chain that exhibits visible effects on 
the intermediate scales. This is due to the presence 
of two full-fledged PS stages. 



D. Yukawa terms 



The effects of the Yukawa couplings can be at lead- 
ing order approximated by constant negative shifts 
of the one-loop coefficients a; — > a[ — a,i + Aa-i 




7.5 8.0 8.5 9.0 9.5 10.0 10.5 11.0 
(a) Chain la 




10.0 10.2 10.4 10.6 10.8 11.0 11.2 11.4 
(b) Chain Va 



n 2 nu 

17 r 




10 11 12 13 

(c) Chain Vila 

FIG. 3: Example of chains with sizeable 126 effects (long- 
dashed curves) on the position of the intermediate scales. The 
solid curves represent the two- loop results in Fig.[l] The most 
dramatic effects appear in the chain la, while moderate scale 
shifts affect chain Va (both "unstable" under small variations 
of the /3-functions). Chain Vila, due to the presence of two 
PS stages, is the only "stable" chain with visible c/i 126 effects. 



with 



Aai 



(4tt) : 



rjjikTrYkY, 



(37) 



The impact of Acn on the position of the unifica- 
tion scale and the value of the unified coupling can 
be simply estimated by considering the running in- 
duced by the Yukawa couplings from a scale t up 
to the unification point [t = 0) . The one-loop result 
for the change of the intersection of the curves corre- 
sponding to a^ 1 (i) and aj(t) reads (at the leading 
order in Aa^): 



Atj 



2tt 



Aai — Aa 
(at - a,) 2 



(38) 
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and 



E. The privilege of being minimal 



Ac 



(a, + aj)(Aaj — Aoj) 



a 3 f 



x [aj\t) 



(39) 



for any i ^ j. For simplicity we have neglected the 
changes in the a,i coefficients due to crossing inter- 
mediate thresholds. It is clear that for a common 
change Aai = Aaj the unification scale is not af- 
fected, while a net effect remains on ajj 1 . In all 
cases, the leading contribution is always propor- 
tional to olJ (t) — (t) (this holds exactly for 
Atu). 

In order to assess quantitatively such effects we 
shall consider first the SM stage that accounts for a 
large part of the running in all realistic chains. The 
case of a low n\ scale leads, as we explain in the fol- 
lowing, to comparably smaller effects. The impact 
of the Yukawa interactions on the gauge RGEs is 
readily estimated assuming only the up-type Yukawa 
contribution to be sizeable and constant, namely 

Tr Y V Y V ~ 1. This yields Aa, 6 x lCT 3 ?/^, 

where the values of the ym coefficients are given 
in Table IXII For i = 1 and j = 2 one obtains 
Aoi ~ -1.1 x 1CT 2 and Aa 2 ~ -0.9 x 10" 2 respec- 
tively. Since af M = § and af M = -f , the first 
term in (f3"9")l dominates and one finds Aa" 1 ~ 0.04. 
For a typical value of ay 1 ~ 40 this translates into 
Acty 1 fcty 1 ~ 0.1%. The impact on tjj is indeed tiny, 
namely Any ~ -1 x 10~ 2 . In both cases the esti- 
mated effect agrees to high accuracy with the actual 
numerical behavior we observe. 

The effects of the Yukawa interactions emerging at 
intermediate scales (or of a non-negligible Tr Yd Yq 
in a two Higgs doublet settings with large tan /3) can 
be analogously accounted for. As a matter of fact, 
in the 5O(10) type of models Tr Y N ~ Tr Y v Y v 
due to the common origin of Yu and Yv. The uni- 
fied structure of the Yukawa sector yields therefore 
homogeneous A<n factors (see the equality of J2k V ik 
in Table IXlf . This provides the observed large sup- 
pression of the Yukawa effects on threshold scales 
and unification compared to typical two-loop gauge 
contributions. 

In summary, the two-loop RGE effects due to 
Yukawa couplings on the magnitude of the unifi- 
cation scale (and intermediate thresholds) and the 
value of the GUT gauge coupling turn out to be 
very small. Typically we observe negative shifts at 
the per- mil level in both nu and ajj, with no rel- 
evant impact on the gauge-mediated proton decay 
rate. 



With all the information at hand we can finally 
approach an assessment of the viability of the vari- 
ous scenarios. As we have argued at length, we can- 
not discard a marginal unification setup without a 
detailed information on the fine threshold structure. 

Obtaining this piece of information involve the 
study of the vacuum of the model, and for 5*0(10) 
GUTs this is in general a most challenging task. In 
this respect supersymmetry helps: the superpoten- 
tial is holomorphic and the couplings in the renor- 
malizable case are limited to at most cubic terms; 
the physical vacuum is constrained by GUT-scale 
F- and £>-flatness and supersymmetry may be ex- 
ploited to studying the fermionic rather than the 
scalar spectra. 

It is not surprising that for non-supersymmetric 
50(10), only a few detailed studies of the Higgs po- 
tential and the related threshold effects (see for in- 
stance Refs. [H, |H, Hi], HH [36[) are available. In 
view of all this and of the intrinsic predictivity re- 
lated to minimality, the relevance of carefully scruti- 
nizing the simplest scenarios is hardly overstressed. 

The most economical 5O(10) Higgs sector in- 
cludes the adjoint 45//, that provides the breaking of 
the GUT symmetry, either 16// or 126//, responsible 
for the subsequent B — L breaking, and 10//, partic- 
ipating to the electroweak symmetry breaking. The 
latter is needed together with 16// or 126// in order 
to obtain realistic patterns for the fermionic masses 
and mixing. Due to the properties of the adjoint rep- 
resentation this scenario exhibits a minimal number 
of parameters in the Higgs potential. In the cur- 
rent notation such a minimal non-supersymmetric 
5O(10) (MSO10) GUT corresponds to the chains 
VIII and XII. 

From this point of view, it is quite intriguing that 
our analysis of the gauge unification constraints im- 
proves the standing of these chains (for Xlla dramat- 
ically) with respect to existing studies. In particu- 
lar, considering the renormalizable setups (126//), 
we find for chain Villa, n\ < 9.1, nu = 16.2 and 
a^ 1 = 45.4 (to be compared to n\ < 7.7 given in 
Ref. @). This is due to the combination of the up- 
dated weak scale data and two loop running effects. 
For chain Xlla we find n\ < 10.8, nu = 14.6 and 
Ojj 1 = 44.1, showing a dramatic (and pathological) 
change from n\ < 5.3 obtained in [9J. Our result sets 
the B — L scale nearby the needed scale for realistic 
light neutrino masses. 

We observe non-negligible two-loop effects for the 
chains VHIb and Xllb (16//) as well. For chain VI- 
Ilb we obtain n\ < 10.5, nu = 16.2 and a^ 1 = 45.6 
(that lifts the B — L scale while preserving njj well 
above the proton decay bound Eq. ([36]) ). A similar 
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shift in n\ is observed in chain Xllb where we find 
rii < 12.5, njj — 14.8 and a v l = 44.3. As we have 
already stressed one should not too readily discard 
njj = 14.8 as being incompatible with the proton de- 
cay bound. We have verified that reasonable GUT 
threshold patterns exist that easily lift nu above the 
experimental bound. For all these chains D-parity is 
broken at the GUT scale thus avoiding any cosmo- 
logical issues (see the discussion in Sect. HI)) . 

As remarked in Sect. lIVBl the limit n\ — n2 leads 
to an effective two-step 50(10) -> G2 -> SM sce- 
nario with a non-minimal set of surviving scalars 
at the G2 stage. As a consequence, the unifica- 
tion setup for the MSO10 can be recovered (with the 
needed minimal fine tuning) by considering the limit 
"2 = nu in those chains among I to VII where Gl is 
either 2 L 2 R 1 X 3 C or 2 L 1 R A C (see Table IJ). From in- 
spection of Figs. [T]-[5] and of Table HTT1 one reads the 
following results: for 50(10) — ► 2 L 2 J? 1 X 3 C -> SM 

45 

we find 

n\ = 9.5, ?i(7 = 16.2 and a^ 1 = 45.5, in case a 
and 

Tii = 10.8, nu = 16.2 and a^ 1 — 45.7, in case b , 
while for 5O(10) — ► 2 L l fl 4 c -> SM 

45 

ri\ = 11.4, nu — 14.4 and a^ 1 = 44.1, in case a 
and 

n\ = 12.6, nu = 14.6 and a^ 1 = 44.3, in case b . 

We observe that the patterns are quite similar 
to those of the non-minimal setups obtained from 
chains VIII and XII in the n% = n% limit. Adding 
the </> 126 multiplet , as required by a realistic matter 
spectrum in case a, does not modify the scalar con- 
tent in the 2l2^1jc 3 C case: only one linear combina- 
tion of the Wh and 126 h bidoublets (see Table Hlf is 
allowed by minimal fine tuning. On the other hand, 
in the liX^c case, the only sizeable effect is a shift 
on the unified coupling constant, namely a^ 1 = 40.7 
(see the discussion in Sect. IIV C| . 

In summary, in view of realistic thresholds effects 
at the GUT (and B — L) scale and of a modest fine 
tuning in the see-saw neutrino mass, we consider 
both scenarios worth of a detailed investigation. 



V. OUTLOOK 

We presented an updated and systematic two- loop 
discussion of non-supersymmetric 5O(10) gauge uni- 
fication with two (and one) intermediate scales. We 
completed and corrected existing analyses by includ- 
ing a thorough discussion of U(l) mixing, which af- 
fects the gauge running already at the one-loop level 



in a number of interesting 5O(10) breaking chains. 
We assessed the relevance of additional Higgs multi- 
plets, needed at some of the intermediate stages in 
order to reproduce a realistic fermionic mass spec- 
trum. Finally, we found and fixed several discrepan- 
cies in the two-loop /3-coefficients. 

The updated results have a non-negligible impact 
on the values of the unification and B — L scales (as 
well as on the value of the unified gauge coupling). 
This is due to the combined effects of the one-loop 
dynamics corresponding to the U(l) gauge mixing 
and of the two-loop RGE contributions. 

We discussed the viability of the different 50(10) 
scenarios on the basis of proton decay and the see- 
saw induced neutrino mass. We were lead to focus 
our attention on the minimal 50(10) setup, emerg- 
ing from a balance of minimal dimensionality Higgs 
representations and a minimal number of parameters 
in the scalar potential. Such a scenario invokes, in 
addition to a complex 10# , one adjoint 45h together 
with one 126# or 16# at the effective level. 

Although the updated values of the unification or 
B — L scales are in some of the setups still conflicting 
with the experimental requirements, they are close 
enough that reasonable spreads in the GUT thresh- 
olds (or a moderate fine tuning in the neutrino mass 
matrix) can easily restore the agreement. This may 
entail the detailed study of the scalar potential of 
the model beyond the tree approximation, that is 
a rather non-trivial task. Nevertheless, the appeal 
of minimality (with supersymmetry confined to the 
Planck scale) motivates us to pursue this study. 
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APPENDIX A: ONE- AND TWO-LOOP BETA 
COEFFICIENTS 

In this appendix we report the one- and two-loop 
/^-coefficients used in the numerical analysis. The 
calculation of the U{X) mixing coefficients and of the 
Yukawa contributions to the gauge coupling renor- 
malization is detailed in Apps. IA II and IA 21 respec- 
tively. 
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TABLE V: The ai and bij coefficients due to pure gauge interactions are reported for the G2 chains with 126^ (left) and 16^ 
(right) respectively. The two-loop contributions induced by Yukawa couplings are given in Appendix IA 21 
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TABLE VI: The and coefficients due to gauge interactions are reported for the Gl chains I to VII with 126^f (left) and 
16 h (right) respectively. The two-loop contributions induced by Yukawa couplings are given in Appendix IA 21 
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TABLE VII: The ai and bij coefficients due to purely gauge interactions for the Gl chains VIII to XII are reported. For 
comparison with previous studies the /3-coefficients are given neglecting systematically one- and two-loops £/(l) mixing effects 
(while all diagonal (7(1) contributions to abelian and non-abelian gauge coupling renormalization are included). The complete 
(and correct) treatment of £/(l) mixing is detailed in Appendix I A II 
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TABLE X: One- and two-loop additional contributions 
to the /3-coefficients related to the presence of the 126 
scalar multiplets in the 2l2_r4 (top) and 2l1_r4 (bottom) 
stages. 
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TABLE IX: The rescaled two-loop /3-coefficients bij com- 
puted in this paper are shown together with the corresponding 
equations in Ref. 0. For the purpose of comparison Yukawa 
contributions are neglected and no U(l) mixing is included in 
chain VIIIa/Gl. Care must be taken of the different order- 
ing between abelian and non-abelian gauge group factors in 
Ref. 0] . We report those cases where disagreement is found 
in some of the entries, while we fully agree with the Eqs. A9, 
All and A16. 
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1. Beta-functions with (7(1) mixing 

The basic building blocks of the one- and two-loop 
/3-functions for the abelian couplings with U(l) mix- 
ing, c.f. Eqs. (|14[) - (Ti"5"|) , can be conveniently written 
as 

9k a gkb = 9sa^sr9rb (Al) 

and 

gkagkbgi c = g sa Ff)grb , (A2) 



where and are functions of the abelian 
charges Ql and, at two loops, also of the gauge cou- 
plings. In the case of interest, i.e. for two abelian 
charges U(1)a and U(1)b, one obtains 

rft = (Qk) 2 , 

^AB — ^BA = QkQk > (A3) 
^BB — (Qk) 2 ' 

and 



^aa = (Qt)\g 2 AA + 9 2 ab) + 2(Qi) 3 Qk(gAAg B A + 9ab9bb) + {Qkf{Qkf{g 2 BA + g 2 BB) , 

r [ A ] B = ^ { b ] a = {QkfQk{g 2 AA + g 2 AB ) + nQi) 2 (Qk) 2 (gAAg B A + gABg BB ) + Qi(Q^) 3 (g 2 B A + 9%b) , (A4) 

r { BB = (Qt) 2 (Qk) 2 (g 2 AA + g\ B ) + 2Qi(Q^) 3 ( 9AA g BA + g AB gBB) + (Qkf(g 2 BA + gls) ■ 



All other contributions in Eq. (fl4|) and Eq. (fT5|) 
can be easily obtained from Eqs. (|A3[) - (|A4[) by in- 
cluding the appropriate group factors. It is worth 
mentioning that for complete £0(10) multiplcts, 

(Qfe ) n (Qk ) m = for n and m odd ( with n + m = 2 
at one-loop and n + m = 4 at two-loop level) . 



By evaluating Eqs. (|A3|) - (|A4|) for the particle con- 
tent relevant to the 2l\ r 1x3c stages in chains VIII- 
XII, and by substituting into Eqs. (fl"4"f - (fl"5")) . one 
finally obtains 



• Chains VIII-XII with 126# in the Higgs sector: 
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• Chains VIII-XII with 16h in the Higgs sector: 
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By setting j X r = Irx = and g X R = gRx = 
in Eqs. (|A5[) - (|A6[) one obtains the one- and two-loop 
/3-coefficients in the diagonal approximation, as re- 
ported in Table [VTTl The latter are used in Figs. [[HI] 
for the only purpose of exhibiting the effect of the 
abelian mixing in the gauge coupling renormaliza- 
tion. 



Yukawa contributions 



The sum can be factorized into the flavour space part 
J2ab Y u b * Y u b = ^[YuY^] times the trace over the 
gauge contractions Tr[AA^] where A = £ £3}. For 
the SM gauge group (with the properly normalized 
hypercharge) one then obtains yiu = j^, y2ii — § 
and ysu = 2, that coincide with the values give n in 
the first column of the matrix (B.5) in Ref. [2lj . 

All of the y p k coefficients as well as the structures 
of the relevant A-tensors are reported in Table IXIl 



The Yukawa couplings enter the gauge /3-functions 
first at the two-loop level, c.f. Eq. © and Eq. (TT4]) . 
Since the notation adopted in Eqs. ©-(O is rather 
concise we shall detail the structure of Eq. ([6]), pay- 
ing particular attention to the calculation of the y p k 
coefficients in Eq. (|2T)]) . 

The trace on the RHS of Eq. ^ is taken over all 
indices of the fields entering the Yukawa interaction 
in Eq. ([7]). Considering for instance the up-quark 
Yukawa sector of the SM the term QlYuUrK + h.c. 
(with h — ia2h*) can be explicitly written as 

YS b e kl 5 z )W Llk U^M + h.c. , (A7) 

where {a, b}, {i,j} and {k,l} label flavour, SU(3) C 
and SU(2)l indices respectively, while 5 n denotes 
the n-dimensional Kronecker S symbol. Thus, the 
Yukawa coupling entering Eq. © is a 6-dimcnsional 
object with the index structure Yjj '£ kl $3j • The con- 
tribution of Eq. (|A7[) to the three y p jj coefficients 
(conveniently separated into two terms correspond- 
ing to the fermionic representations Ql and Ur) can 
then be written as 

y P u = -^[c$\Ql) + c$\u r ) 

x jr Y$ h e kl 8 5 )Y$ b *e kl 8£ (A8) 

ab,ij,kl 



20 



Go 



y P k 



k Gauge structure Higgs representation Tensor A Tr[AA^] 



1y 
2l 
3 C 




U 
D 
E 



QL kj Uhhi 
Ql_ kj D' R h l 
LLkE % B h l 



h> : (+1,2,1) 



^4 

021 



2 L 
3c 



3 

2 
3 
2 

1 fl 

2 V 2 

1 
2 
I 

2 



1 
2 
1 
2 

1 /¥ 

2 y 2 
1/31/3 

" 2 Y 2 2 Y 2 
3 3 
2 2 



u 

D 

N 
E 



QL_ kj D R h l 
L L kN R hi 
LLkEnh 1 



h> : (2, +|,0,1) 



2l 


f 3 


1^ 




2b 


3 


1 


Q 


lx 


1 


3 


L 


3 C 


u 






2 L 






F u 


lx 


2 


:; 


F D 


4c 


U 







Ql k Q c L T<P ln 

j^k j^c miln 



FLkjF R Jl hi 
FLkjFR l h l 



<t> ln : (2,2,0,1) 



h l : (2, +|,1) 



^kl^mn 



12 
4 



2l 
2b 

4c 



-pifc rpc m j Zn 

'l F Lj (/> 



16 



2t 
lx 
4c 

2i 
2 B 
4c 




FL kj Fg l H? 
~K k] F° i H la 



L Lj 



H la : (2, +±,15) 



(2,2,15) (Ta)l 



15 
15 
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TABLE XI: The two-loop Yukawa contributions to the gauge sector /3-functions in Eq. ([20} are detailed. The index p 
in y pk labels the gauge groups while k refers to flavour. In addition to the Higgs bi-doublet from the 10-dimensional 
representation (whose components are denoted according to the relevant gauge symmetry by h and <j>) extra bi- 
doublet components in 126// (denoted by H and <&) survives from unification down to the Pati-Salam breaking scale 
as required by a realistic SM fermionic spectrum. The T a factors are the generators of SU(4)c in the standard 
normalization. As a consequence of minimal fine tuning, only one linear combination of IQh and 12&h doublets 
survives below the SU(4)c scale. The U(1)r,x mixing in the case 2_l1.r1x3c is explicitly displayed. 
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